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QJ , Abstract 

^ ' We present a measurement of the inelastic, non diffractive J ftp photoproduction 

cross section in the reaction e^p — > e^J/ipX with the ZEUS detector at HERA. 
The J/'0 was identified using both the /x~''/U~ and e+e^ decay channels and events 
^ ■ were selected within the range 0.4 < z < 0.9 (0.5 < z < 0.9) for the muon (electron) 

decay mode, where z is the fraction of the photon energy carried by the J/ip in 
the proton rest frame. The cross section, the and the z distributions, after 
having subtracted the contributions from resolved photon and diffractive proton 
dissociative processes, are given for the photon-proton centre of mass energy range 
50 < < 180 GeV; pj^ is the square of the 3/ip transverse momentum with respect 
to the incoming proton beam direction. In the kinematic range 0.4 < z < 0.9 and 
> 1 GeV^, NLO calculations of the photon-gluon fusion process based on the 
colour-singlet model are in good agreement with the data. The predictions of a 
specific leading order colour-octet model, as formulated to describe the CDF data 
on J/ip hadroproduction, are not consistent with the data. 
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1 Introduction 



The inelastic reaction e~^p —>■ e'^J/tpX in the photoproduction regime [Q^ ~ GeV^, 
where is the photon virtuahty) is thought to proceed via direct photon-gluon fusion, 
diffractive proton dissociation or resolved photon processes. These three possibilities are 
shown in Fig. [I|. In this paper we are primarily interested in the contribution from the 
direct photon-gluon fusion process, shown in Fig. la, for which full next-to-leading order 
(NLO) QCD calculations are available in the framework of the colour-singlet model 0. 
The predicted cross section is sensitive to the gluon density in the proton at an energy scale 
corresponding approximately to the heavy quark mass. Perturbative QCD calculations 
using as input the parton densities extracted from other processes can therefore provide 
a consistency test for QCD. In inelastic J/ip photoproduction the concepts of direct and 
resolved photon contributions remain distinct up to the NLO level. This is due to the 
particular spin and colour state projection involved in the calculation, which makes this 
reaction different from those involving light and open heavy quark production, in which 
only the sum of direct and resolved terms is unambiguously defined at NLO . 

The different processes in Fig. |I| can be distinguished by means of the inelasticity variable 
z, the fraction of the photon energy carried by the J/ip in the proton rest frame @]. The 
diffractive proton dissociation process dominates the high z region {z > 0.9), the photon- 
gluon fusion the intermediate z region, while the resolved photon process is dominant 
in the low z region. There are also differences in the pt distribution, pt being the J/ip 
transverse momentum with respect to the incoming proton beam direction. Compared to 
photon-gluon fusion and resolved photon processes the diffractive reaction produces J/ip 
mesons with relatively low px {pr ^ 1 GeV) . 

The direct photon-gluon fusion process is described in in the framework of the colour- 
singlet model [0], in which the initial state photon and gluon interact giving a final state 
cc pair with the J/ip quantum numbers through the emission of a hard gluon in the final 
state (7 + — > J/V' + 92)- When a similar model was used to study J/ip hadroproduction 
at the Tevatron, the predictions, at lowest order in the strong coupling constant a^, 
underestimated the data by about one order of magnitude. The measured cross section 
could be explained in part by adding colour-octet contributions 0. In this case the cc 
pair is produced in a colour-octet state (short distance process) and later binds to form 
Si J/ijj (long distance process). While the short distance terms are calculable through 
perturbative QCD, the long distance terms are nonperturbative and have to be determined 
from the data themselves. It is therefore interesting to look for evidence of the octet 
mechanism at HERA, where it is expected to contribute at high z H. 



Previous fixed target experiments both in the photoproduction and in the electro- 



production [jTT], |T^, |T4[ regime measured the inelastic cross section for photon- 



proton centre of mass energies, W, between 10 and 20 GeV. The HI Collaboration [|T5[ 
has published results on inelastic J/ip production in the interval 30 < < 150 GeV. 

In the following sections, after a brief description of the experimental conditions, we 
discuss the kinematics of inelastic J/ip production and the criteria used to select events 
in the region where the direct photon-gluon fusion process is dominant. We then evaluate 
the cross section for this process in the range 50 < 1^ < 180 GeV and 0.4 < z < 0.9 
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for the muon case and in 90 < < 180 GeV, 0.5 < z < 0.9 for the electron case. The 
cross section is extrapolated to ^ = assuming the direct photon-gluon fusion model. 
Comparisons with NLO calculations are discussed in section ^ for the restricted kinematic 
range z < 0.8 and p\ > 1 GeV^, where the calculations are reliable. 

The data were collected in 1994 and correspond to an integrated luminosity of 2.99±0.05 
pb-i. 



2 Experimental Conditions 
2.1 HERA 

During 1994 HERA operated with a proton beam energy of 820 GeV and a positron 
beam energy of 27.5 GeV. There were 153 colliding proton and positron bunches together 
with an additional 17 unpaired proton bunches and 15 unpaired positron bunches. The 
root mean square (rms) proton bunch length was approximately 20 cm while the positron 
bunch length was small in comparison. The time between bunch crossings was 96 ns. The 
typical instantaneous luminosity was 1.5 x 10^° cm~^ s~^. 



2.2 The ZEUS Detector 

The main ZEUS detector components used in this analysis are outlined below. A detailed 
description of the ZEUS detector can be found elsewhere [|16|. In the following the ZEUS 



coordinate system is used, the Z axis of which is coincident with the nominal proton 
beam axis, the X axis is horizontal and points towards the centre of HERA and the Y 
axis completes a right-handed coordinate system. The origin of the coordinates is at the 
nominal interaction point. 

The momenta and trajectories of charged particles are reconstructed using the vertex 
detector (VXD) Q and the central tracking detector (CTD) The VXD and the 

CTD are cylindrical drift chambers covering the angular region 15° < 6 < 164° (where 6 
is the polar angle with respect to the proton direction). The chambers are located in a 
magnetic field of 1.43 T produced by a thin superconducting solenoid. 



The high resolution uranium-scintillator calorimeter (GAL) |TP[ surrounding the coil is 
divided into three parts, the forward calorimeter (FCAL), the barrel calorimeter (BGAL) 
and the rear calorimeter (RGAL). They cover polar angles from 2.6° to 36.7°, 36.7° to 
129.1°, and 129.1° to 176.2°, respectively. Each part consists of towers which are longitu- 
dinally subdivided into electromagnetic (EMG) and hadronic (HAG) readout cells. The 
GAL also provides a time resolution of better than 1 ns for energy deposits greater than 
4.5 GeV, and this timing is used for background rejection. 

The hadron electron separator (HES) consists of silicon detectors 400 fim thick. In 
the 1994 running period only the rear part (RHES) was operational. The RHES is located 
in the RGAL at a depth of 3.3 radiation lengths, covering an area of about 10 m^. Each 
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silicon pad has an area of 28.9 x 30.5 mm^, providing a spatial resolution of about 9 mm 
for a single hit pad. If more than one adjacent pad is hit by a shower, a cluster consisting 
of at most 3x3 pads around the most energetic pad is considered. This allows a more 
precise reconstruction of the position with a resolution of about 5 mm for energies greater 
than 5 GeV. The RHES measures the energy deposited by charged particles near the 
maximum of an electromagnetic shower. 



The muon detectors |21| , situated outside the calorimeter, consist of limited streamer tube 
(LST) planes with the inner chambers in front of the magnetised iron yoke and the outer 
chambers behind it. Owing to the low momentum of the J/V' decay muons, only the 
inner chambers (BMUI and RMUI) were used in the present analysis. The BMUI and the 
RMUI cover the polar angular ranges 34° < 9 < 135° and 134° < 9 < 171°, respectively. 

The luminosity is determined from the rate of events due to the Bethe-Heitler process 
e+p — > e"*"7p, where the photon is measured by the calorimeter of the luminosity detector 



(LUMI) located in the HERA tunnel in the direction of the outgoing positron beam p2 
For the measurements presented in this paper the luminosity was determined with a 
precision of 1.5%. 



3 Kinematics 



Schematic diagrams for the reaction: 

where each symbol in parentheses denotes the four-momentum of the corresponding par- 
ticle, are shown in Fig. |I|. 

The kinematics of the inclusive scattering of unpolarised positrons and protons is described 
by the positron-proton centre of mass energy squared (s) and any two of the following 
variables: 



• Q"^ = —q^ = —{k — k'Y, the negative four-momentum squared of the exchanged 
photon; 

• y = {q ' ■ P), the fraction of the positron energy transferred to the hadronic 
final state in the rest frame of the initial state proton; 

• W"^ = + = + 2y(/c ■ P) + ^ ys, the centre of mass energy squared of 
the photon-proton system, where Mp is the proton mass. 



Restricting our measurement to photoproduction events where the outgoing positron is 
not in the CAL acceptance, the value ranges from the kinematic minimum Q^j^ = 
Mgt/^/(l — y) ~ 10^^° GeV^, where Me is the electron mass, to the value at which the 
scattered positron starts to be observed in the uranium calorimeter, Q^ax ~ 4 GeV^. The 
median is approximately 10~^ GeV^. 
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The value of y was determined by the Jacquet-Blondel formula 



= (2) 



-'e 



where the sum runs over all the calorimeter cells and Ei is the cell energy, is equal 
to Ei cos 6i, where 6i is the polar angle of the cell measured with respect to the nominal 
vertex, and Ee is the incoming positron energy. The value of W is determined from the 
relation W"^ = VjbS- 

For reaction (|ID the inelasticity variable, is defined by 

_ P ■ pj/y, _ {Ej/^-pzj/^) 
P-q - 2yE, ' 

where Ej/^ is the J/ip energy and PZj/^ is its momentum component along the Z direction. 

In the proton rest frame, z is equal to where E^ is the photon energy. Experimentally 
z is estimated from: 

2yjBE, ~ Ei{E^~pzJ ■ ^ ^ 

In the estimations of yjB and z, using equations (2) and (4), the contribution of the two 
leptons from the J/ip decay was accounted for by including in the sum their momenta as 
measured in the central tracking detectors while discarding their calorimetric deposits. 



4 Event Selection 



The selection of the muon and electron decay channels followed different paths, except 
for the common veto requirements at the first level trigger (FLT), which reject proton- 
gas background events occuring upstream of the nominal interaction point and which are 
therefore out of time with respect to the e~^p interactions. 



4.1 Muon Mode 

The candidates for the J/ip ^ f^^f^~ channel were selected using the three level ZEUS 
trigger system. At the FLT a coincidence between track segments in the CTD, energy 
deposits in the CAL and hits in the BMUI or RMUI was used to select muon candidates. 
The CAL was divided in Z — (p regions (0 being the azimuthal angle around the Z axis) 
associated with the corresponding zones of the barrel and rear muon chambers. A signal 
above threshold in one of the CAL regions in conjunction with a hit in the associated barrel 
or rear muon chamber defined a CAL-BMUI/RMUI match. This regional matching was 
demanded together with the requirement of tracks in the CTD pointing to the nominal 
vertex. 

At the second level trigger (SLT), the total energy in the calorimeter {Etoi = ^iEi) and 
the Z component of the momentum (Sjp^. = Sji?jCos6'j) were calculated. The sums 
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run over all calorimeter cells i with an energy, Ei, above threshold at a polar angle, 9i, 
measured with respect to the nominal vertex. In order to remove proton-gas interactions, 
events with the ratio T,ipzJ Etoi greater than 0.96 were rejected. Part of the cosmic ray 
background was rejected at the SLT by using the time difference of the energy deposits 
in the upper and the lower halves of the calorimeter. 

At the third level trigger (TLT) a muon candidate was selected when a track found in 
the CTD matched both a cluster with a calorimeter energy deposit consistent with the 
passage of a minimum ionising particle (a m.i.p. cluster 0) and a track in the inner muon 
chambers. An event containing a muon candidate in the rear (barrel) region was accepted 
if the (transverse) momentum of the CTD track exceeded 1 GeV. 

The TLT algorithm was again applied in the offline analysis, but now the results of the 
full event reconstruction were used. The tracks corresponding to the two muons from 
the J/ip decay had to satisfy the following criteria, where the subscript 1 denotes the 
triggering muon and the subscript 2 the other muon and p indicates the momentum of a 
muon and pt its transverse momentum: 

• Pi > 1 GeV (rear region); pn > 1 GeV (barrel region); 

• P2>1 GeV; 

• Pti + Pt2 > 2.8 GeV; 

• pseudorapidities0 |?7i_2| < 1.75; 

• the second muon track has to match a m.i.p. cluster in the GAL. 

Gosmic rays were rejected by requiring that the two muon tracks were not coUinear: 
events with Q > 174° were rejected, where Q is the angle between the two tracks at the 
interaction point. 

The final inelastic data sample was defined by requiring an energy deposit greater than 
1 GeV in a cone of 35° around the forward direction (excluding the calorimeter deposits 
due to the muons). Elastically produced J/ip mesons were thus excluded. The data were 
further restricted to the W interval 50 to 180 GeV where the acceptance is above 10%. 
The data sample was divided into the three categories: 

• events with z in the interval 0.9 to 1; 

• events with z in the interval 0.4 to 0.9; 

• events with z < 0.4. 

The first category is interpreted as coming mainly from the diffractive proton dissociation 
process. The second one is dominated by the photon-gluon process (direct process) and 
the third is a combination of direct and resolved processes. The fi~^^~ invariant mass for 

cluster is defined as a group of contiguous cells in the CAL with energy above a set threshold. 
^The pseudorapidity is defined as 77 = — lntan(|). 
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the second category is shown in Fig. ga, fitted with a Gaussian plus a fiat background 
giving a mass of 3.086±0.004 GeV. The rms width is 39±4 MeV, consistent with the 
Monte Carlo expectations. The invariant mass distribution for z < OA events is shown in 
Fig. Table ^ contains the fitted number of events above background for each category 
and for various W ranges. The data so collected correspond to events with < 4 GeV^. 
The events selected in the chosen W range have J2i{Ei —pzj < 20 GeV. Events with the 
scattered positron in the GAL are expected to have a J2i{Ei — pzj ~ 2Ee = 55 GeV. A 
cross check with an electron finder confirmed the absence of large events in the sample. 

4.2 Electron Mode 

Inelastic J/ip ^ e^e~ candidates were triggered at the FLT by demanding the two con- 
ditions: 

1. at least one of the following requirements on the GAL energies: 

• GAL total energy > 15 GeV; 

• GAL-EMG energy > 10 GeV; 

• GAL total transverse energy > 11 GeV; 

• BGAL-EMG energy > 3.4 GeV; 

• RGAL-EMG energy > 2.0 GeV; 

2. at least one GTD track associated with the nominal vertex. 

At the SLT events were rejected if J^iPZi/Exot was greater than 0.92 (with and ETot 
defined as in the previous section). In addition only events satisfying the following two 
conditions were accepted: 

• — pzi) > 4 GeV, where the sum runs over all calorimeter cells i; 

• the sum of the total energy deposits in BGAL-EMG and RGAL-EMG was greater 
than 3 GeV. 

At the TLT a fast electron identification was carried out by using information from the 
GTD and GAL. Clusters were identified as electrons if at least 90% of the cluster's energy 
was deposited in the electromagnetic section. The tracks from the GTD were extrapolated 
towards the GAL and matched to the nearest cluster within 30 cm of the extrapolated 
track at the GAL face. An event was accepted if at least two oppositely charged tracks, 
identified as electrons, were found each with a momentum exceeding 0.5 GeV and a 
transverse momentum greater than 0.4 GeV; in addition, the two tracks were required to 
originate from points less than 7 cm apart along the Z axis. The invariant mass of the 
track pair, assuming the electron mass for each track, had to be greater than 2 GeV. 
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The initial offline selection was based on the TLT track-cluster matching algorithm, but 
using the full tracking and CAL information. 

Since the transverse momentum {pt) spectrum of the background tracks peaks at low pt 
values, both electron track candidates were required to have pt greater than 0.8 GeV. Also, 
both tracks had to originate from the event vertex and satisfy the condition \ri\ < 1.75. 
The large background coming mainly from low energy pions faking electrons was further 
reduced firstly by requiring a tighter matching of the tracks to the electromagnetic clusters, 
with a track-cluster separation at the CAL face less than 25 cm, secondly by demanding 
clusters with small longitudinal and radial dimensions and thirdly by imposing a cut 0.4 
< Eciuster/p < 1-6, where Eduster is the energy of the electromagnetic cluster and p is 
the momentum of the associated track. This cut was chosen since, for electrons in the 
momentum range 1-3 GeV (typical for electrons from J/ip decays in the present analysis), 
the inactive material in front of the calorimeter means that the Eciuster/p ratio is about 
0.8 with 20% resolution. 

A significant reduction in the remaining background was achieved by using the RHES 
and the information on the specific ionisation energy loss, dE/dx, as evaluated from the 
CTD. The dE/dx of a track, calculated from the truncated mean of a distribution of pulse 
amplitudes where the lowest 30% and the highest 10% were discarded, had a resolution 
of about 12%, averaged over a broad range in t] {\r]\ < 1.5). By requiring dE/dx for one 
of the electron candidate tracks to be consistent with that expected for an electron, 93% 
of the e+e~ pairs were retained, while discarding two thirds of the background^. Because 
the identification of electrons via dE/dx is not well understood for low angle tracks, the 
use of dE/dx was limited to tracks with Ir^l < 1. For rj < —1 the RHES can be used for 
electron identification. A HES electron cluster was defined as a group of adjacent silicon 
pads each with an energy deposit above 0.6 m.i.p., and the total energy of all pads above 
5 m.i.p. Track-HES cluster matching was then performed for tracks already matched to 
a CAL cluster, requiring that the distance between the HES cluster and the extrapolated 
track be less than 10 cm. The efficiency of the cut on the RHES cluster energy above 
5 m.i.p. was estimated to be 75% using an almost background-free J/ip — > e"'"e~elastic 
sample. 

Events were then accepted if the electron tracks satisfied one of the following requirements: 
if both tracks lie in the range \ri\ < 1 they had to satisfy the aforementioned dE/dx cut; if 
one track was in the range \t]\ < 1 and the other in r] < —1 they had to satisfy the dE/dx 
and the RHES cuts, respectively. All the other track combinations were not considered 
due to the presence of high background. 

The final inelastic data sample was defined by requiring an energy deposit greater than 
1 GeV in a cone of 35° around the forward direction (excluding the calorimeter deposits 
due to the electrons). A minimum value of z was required {z > 0.5) to avoid the low z 
region which is dominated by large background. The data were also restricted to the W 
range 90 to 180 GeV, where the acceptance is high. Using the z variable the electron data 
sample was divided in two categories: z > 0.9 and 0.5 < z < 0.9. 

^The 93% efficiency was computed from 7 conversions and almost background-free J/ip — s- e+e~elastic 
events. 
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Figure ^ shows the mass distribution of the electron pairs for the second category. A 
clear peak at the J/ip mass is observed. The solid line shows an unbinned likelihood fit 
in which a Gaussian resolution function has been convoluted with a radiative J/ip mass 
spectrum and a linear distribution to describe the background (dashed line). The mass 
estimated by the fit is 3.089 ± 0.010 GeV. The rms width is 40 ± 9 MeV, consistent 
with the MC expectation. Table |I] contains, for the two categories and W ranges, the 
fitted number of events above background. As for the muon sample, the data so collected 
correspond to events with < 4 GeV^. 



5 Monte Carlo Simulation and Acceptance Calcula- 
tion 

Inelastic J/ip production from direct photon-gluon fusion was simulated using the colour- 



singlet model as implemented in the HERWIG |2J] parton shower generator. The range 
of was from the kinematic limit (~ 10"^'' GeV^) to 4 GeV^. The energy scale, /i^, at 
which the gluon distribution is evaluated was chosen to be /i^ = 2stu/ {s^ + 1"^ + u'^), where 
s, t and u are the Mandelstam variables of the photon-gluon fusion process. The mean 
value of fi^ is 7 GeV^. The gluon structure function in the proton was parameterized with 
the MRSD'_ p5i distribution. 



For resolved 3/ip production the PYTHIA parton shower generator was used with 



the GRV proton |^ and photon parton densities. The matrix elements for resolved 
photon processes were computed in the colour-singlet framework. 

Production of J/ip mesons accompanied by diffractive proton dissociation was simulated 
with EPSOFT |]2D|. This generator is based on the assumption that the diffractive cross 
section is of the form da / d\t\dMff oc e~^'''*'/M^, where is the mass of the dissociative 
system, and t is the four-momentum transfer squared at the proton vertex. The value of 
bd was chosen to be 1 GeV~^ to reproduce the observed p"^ distribution of events with 
z > 0.9. For the Mn distribution, the value (3 = 2 was used. The simulation of the 
dissociative system includes a parametrisation of the resonance spectrum. 

In the muon case a mixture of HERWIG (78ig)% and EPSOFT events gives the best 
description of the z distribution with z ranging from 0.4 to 1. For the electron case 
the percentage is (79±6)%. The same mixture also describes well all the reconstructed 
kinematic variables (see Fig. ^ for some examples). The resolution in z is 2% ai z = 0.9 
and increases to 12% at z = 0.4. The measured values of z suffer from a systematic 
shift due to the energy loss in the inactive material in front of the calorimeter and to the 
undetected particles escaping in the beampipe. This shift is 20% at z = 0.4 and becomes 
negligible at z = 0.9. The shift of z was corrected using the HERWIG Monte Carlo. 

The acceptance was estimated as the ratio of the number of accepted Monte Carlo photon- 
gluon fusion events to the number generated in the selected kinematic range. The accep- 
tance, calculated in this manner, accounts for the geometric acceptance, for the detector, 
trigger and reconstruction efficiencies, and for the detector resolution. Table |^ reports 
the acceptances in various W ranges determined for each decay mode. 
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6 Backgrounds to the photon-gluon fusion process 



In this section we discuss all the resonant processes which are backgrounds to the photon- 
gluon fusion process. To calculate integrated and differential cross sections the analysis 
was restricted to the region 0.4 < z < 0.9 and 50 < < 180 GeV for the muon mode, and 
to the region 0.5 < z < 0.9 and 90 < < 180 GeV for the electron mode. The upper z 
cut is necessary to exclude diffractive i /ip production; the lower cut restricts the data to a 
range with low background and where the photon-gluon process is expected to dominate 
over the resolved production. Table |I] reports the numbers of events coming from the 
fit to the ^/ip mass peak, divided into four W intervals for the muon mode and three 
intervals for the electron mode. These numbers include proton diffractive dissociation 
events which migrated from the region above z = 0.9. The background from diffractive 
events was estimated to be f'^'^-^-l' = (8 ± 2)% in the muon mode and f'^'^f-f = (4 ± 1)% 
in the electron mode using EPSOFT and following the method explained in section ^. 
The cross sections were corrected for the estimated fraction of proton dissociative events. 
The difference in the size of the contamination between muon and electron decay mode 
reflects the different W range covered. 

At the lower end of the z range one has to consider resolved J ftp photoproduction events 
for z > 0.4(0.5) as well as migration of resolved events from z < 0.4 (0.5) into the region 
studied. In Fig. the invariant mass plot of the muon decay mode for events with z < 0.4 
and 50 < W < 180 GeV is shown. The estimated number of J/ip events is 19 ± 6. As 
determined from the z shape of resolved and direct photon Monte Carlo generated events, 
about 50% of the detected events can be attributed to resolved photoproduction. Their 
contribution in the 0.4 < z < 0.9 range, after diffractive proton dissociation subtraction, 
is 3^2%- This background has been subtracted from the signal events assuming the W 
and z dependence given by PYTHIA Monte Carlo. Due to the small statistics of the 
electron sample the muon result was also used in the electron case. 

Production of ip' mesons with subsequent decay into J/ip is a contribution not included 
in the simulation. It is estimated in []1|, through phase space considerations, to be 15% 
of the J/ip integrated cross section. This result is in good agreement with estimations 



made by using the value of the ^p' to 3 /ip ratio coming from low energy data ||^, |T0[. This 
contribution was not subtracted. 



7 Systematic Errors 

Several factors contribute to the systematic errors in the inelastic J/ip cross section mea- 
surement. In the following they are divided in two categories: decay channel specific errors 
and common systematic errors. The first category contains systematic errors specific to 
the electron or muon decay channel, while the second contains systematic errors common 
to both decay modes. 
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Decay channel specific errors: 



• Trigger, for the muon mode the principal source of uncertainty is the FLT calorime- 
ter trigger. The corresponding error was estimated by using independent muon 
triggers, which use different calorimeter trigger logic or do not use the calorime- 
ter at all. Also for the electron mode the uncertainty is dominated by the FLT 
calorimeter simulation and the corresponding error was estimated by using inde- 
pendent calorimeter triggers. The size of the error depends on the W range and is 
of the order of ±4% in the lowest W bin and of ±1% in the highest one. 

• Event selection: this class comprises the systematic errors due to the uncertainties 
in the measurement of momentum, transverse momentum and 7] of the leptonic 
tracks. For the muon channel this class contains also the uncertainties coming from 
the pti + Pt2 > 2.8 GeV cut and the collinearity cut. For the electron channel 
uncertainties in the definition of an electron cluster also contribute. Each cut was 
varied within a range determined by the resolution of the variable to which the cut is 
applied and the different contributions obtained were summed in quadrature. This 
error amounts to ±2%, almost independent of W. 

• Muon chamber efficiency: the systematic error (±2%) attributed to uncertainties in 
the muon chamber reconstruction efficiency was estimated using cosmic ray events. 

• dE/dx: the error takes into account the uncertainty of the variation in the efficiency 
of the dE/dx cut as a function of the track's polar angle, not reproduced by the 
Monte Carlo simulation. The size of the error depends on the W range and is of 
the order of ±7% in the lowest W bin and of ±5% in the highest one. 

• RHES: the error (+1%) was estimated by rising the cut value from 5 m.i.p. to 6 
m.i.p. in the Monte Carlo simulation only. This was done to take into account 
possible differences of the calibration in the simulation and in the data. 

• Fitting procedure: different fitting procedures for the ^/ip mass peak were applied to 
the electron channel and the results were the same within the statistical uncertainty. 
The size of the error is -3%. 

• Branching Ratio: the error on the branching ratio J/ip ^ is used as quoted in 
H (±3.2%). 

Common systematic errors: 



Parton density: the uncertainty in the acceptance resulting from the uncertainty 
in the form of the gluon distribution was estimated by changing the default gluon 
distribution (MRSD'_ [||]) with others (GRV MRSG H, MRSA' Q). This 



error is of the order of ±1' 
it contributes -8%. 



in all the W ranges, except in the highest W bin where 



Energy scale: the calorimeter energy scale in the Monte Carlo was varied by ±5%. 
This affects mainly the W and z determinations and gives a ±2% contribution. 
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Proton dissociation: The parameters [3 and ba were varied from /? = 2 to /? = 2.5 
{p = 2.20 ± 0.03 is the result of ^) and from bd = 0.9 GeV-^ to ba = 1.3 GeV'^ 



(from the analysis of the distribution of the events with z > 0.9). This systematic 
uncertainty is dominated by the number of events with z > 0.9 and is concentrated 
in the range 50 < < 90 GeV, where it gives a contribution of ±5%, while in the 
other bins it is of the order of the ±1% or lower. 

Resolved photon: this systematic error contains contributions from the limited statis- 
tics of the muon sample for z < 0.4 and the uncertainties in the Monte Carlo mod- 
elling of the resolved process. The corresponding error was evaluated to be ±4%, 
independent from W. 

z extrapolation: the cross sections were measured down to z = 0.4 for the muon 
channel and to z = 0.5 for the electron channel and then extrapolated to zero using 
the HERWIG Monte Carlo. The uncertainty (±3.5%) on the extrapolation was 
evaluated by varying Aqcd and the charm mass, rric, in the NLO calculation. 

Angular distribution: The angular distribution of the decay leptons was modelled 
using the form (1 + acos^ 6*), where 6* is the decay angle of the leptons in the J/ip 
rest frame with respect to the direction of the J/ip momentum in the laboratory. 
The data are best described with « = 0, that is a flat distribution. As a systematic 
check a was varied by one standard deviation (i.e. up to a = 0.5); this gave an 
error growing from +5% in the lowest W bin to +8% in the highest W bin. 



• Luminosity: as indicated in section |2.2| , the uncertainty on the luminosity determi- 
nation is ±1.5%. 

The total systematic error, given by the sum in quadrature of all the common and un- 
common systematic errors, is of the order of ±10% (±15%) for the muon (electron) decay 
channel. 



8 Results 



8.1 Cross Section Calculation 

The electroproduction cross section for inelastic J/ip production, after subtracting the 
contributions of diffractive proton dissociation and of resolved photon processes, is calcu- 
lated as: 

aie^p - e+J/'0X) = (5) 

where N^vt denotes the background subtracted number of J/ip signal events, A the accep- 
tance, C the integrated luminosity and B the J/ip leptonic branching fraction namely 



(6.01 ± 0.19)% for /i"*"/! and (6.02 ± 0.19)% for e+e . The photoproduction cross section 
is related to the ep cross section by ||33| 

o^y^jl^x ^- ^- , (6) 
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where a^p^j/^x is the mean cross section in the measured range of W (corresponding to 
the hmits |/mm, Umax) and . The effective flux, of virtual photons from the positron 
is computed as: 

$T = / I Hy,Q')dydQ' = I —--[1 + {1-yf - —^]dydQ\ 



where a is the electromagnetic couphng constant. The integrals run from Q^j„ = M^-^ 

to Qliax = 4 GeV^ and from ymin = W^in/s to ymax = W^ax/^ where Wmin and Wmax 
are the minimum and maximum values of W, respectively, in each chosen interval. 

The electro- and photoproduction cross sections are summarized in Tab. ^ for the two 
J/ip decay channels. The first error is statistical and the second comes from adding in 
quadrature all the systematic errors described in section |^. The cross sections measured 
in the restricted z ranges were extrapolated to 2; = (using HERWIG) in order to be 
able to compare with other available data. The size of the extrapolation (~ 10% for the 
muon decay mode and ~ 25% for the electron decay mode) and the associated systematic 
error are shown in Tab. ^ and given in section |^. The difference in the size of the 
extrapolation between muon and electron decay modes is due to the different z ranges 
measured. Table |^ reports for the lowest W bin the photoproduction cross section from 
the muon channel only and for the following three W bins the combined muon and electron 
photoproduction cross sections. These results are given both in the range 0.4 < z < 0.9 
and in the extrapolated range z < 0.9. To obtain the combined results a weighted mean 
was calculated; the weights were obtained by summing the statistical and decay channel 
specific errors in quadrature. The first error for the combined results in Tab. ^ is the error 
on the weighted mean, the second is given by the sum of the common systematic errors 
added in quadrature. The photoproduction cross sections for z < 0.9 are shown in Fig. ^ 
together with those found by the HI collaboration |jl5[ and a compilation of fixed target 



results |T^, |T^. The ZEUS and HI data are compatible. The cross section rises as W 
increases. The curves in the plot correspond to the NLO calculation |l|] computed with 
no cut on the pt of the J/ip , z < 0.9, with a charm mass (mc) of 1.4 GeV, Aqcd = 300 
MeV and with renormalization and factorization scales of ^J^m^. Dashed, continuous 
and dashed-dotted curves are obtained with different gluon distributions (MRSG [pO| , 



GRV and CTEQ3M respectively), compatible with those extracted from F2 
measurements at HERA |p5| , The theoretical predictions were multiplied by a factor 
1.15 to take into account t/^' production. The predictions are in qualitative agreement 
with the data, but the distinction among different parton densities is not possible because 
the NLO calculation is not well behaved in the limit 0. Furthermore, there is a 

significant dependence on the values of the charm mass and I^qcd- The cross section 
varies as 1/m^ and as illustrated by the dotted line which is calculated with vrtc = 1.55 
GeV and Aqcd = 215 MeV and using the GRV gluon distribution. 

A more quantitative comparison between data and theory can be made in the restricted 
kinematic range p\. > 1 GeV^, where the calculation is much more reliable. The NLO 
computation now allows an absolute comparison between data and models. The cross 
sections for this kinematic range are summarized in Tab. |]for the two ^/ip decay channels 
and in Tab. 1^ for the combined result, with the additional requirement z < 0.8 in order 



to compare with |T5[. The measured cross sections for z < 0.8 and > 1 GeV are 



12 



displayed in Figure ^ The curves represent the NLO calculation using the different gluon 
distributions cited above. Data and theory are in good agreement using rric = lA GeV, 
Kqcd = 300 MeV and v^mc as renormalization factor. While for z < 0.9 and all the 
predictions are significantly different for different parametrizations of the gluon density, 
this is not so in the more restricted domain z < 0.8 and > 1 GeV^. This is a 
consequence of the cut: the gluon distribution is probed at larger Xg values, where the 
differences between the various gluon densities are smaller; here Xg is the proton energy 
fraction carried by the incoming gluon. In the present analysis, with z < 0.9 we explore 
the range 4 ■ 10"^^ Xg ^10^^ and the range 10 Xg ^10 ^ in the restricted interval. 



8.2 Transverse Momentum Distribution 

Figure H shows the differential cross section da/dp'^ for z < 0.9 and 50 < 1^ < 180 
GeV using only the muon sample. The background contributions listed in section ^ were 
subtracted bin by bin. The curve shows the NLO prediction obtained with rric = 1.4 GeV, 
Aqcd = 300 MeV, \/2mc as scale and GRV for the gluon distribution^. For Px > ^ GeV^ 
data and theoretical calculation are in good agreement. A fit of the function 



da 



2 -Ae-'PT (8) 



dprp 

to our data was performed in the range 1 < < 9 GeV^ giving 

b = 0.32 ± 0.03 (stat) ± O.Ol(syst) GeV~^ (9) 

A similar fit to the NLO calculation |l[] yields a slope of b = 0.3 GeV^^ above > 1 
GeV^. The systematic error contains contributions from all the classes of systematic errors 
discussed in section |^ and also from the change in the p"^ fitting interval. 



8.3 Distribution of z 

Figure ^ shows the differential cross section da/dz for p'^ > 1 GeV^ and 50 < < 180 
GeV as obtained using only the muon sample. It is compared to the NLO calculation [|l| 
discussed in section ^]l] and with the parameters used in |8]^. Agreement in shape and 
normalization is found within the errors. Our data are in good agreement with the result 



of the HI collaboration |T5|. 



Recently there has been theoretical activity attempting to solve the discrepancy between 
the J/tp production cross section measurements in hadronic reactions and the colour- 



singlet model by invoking additional octet contributions |^. A specific leading order 
calculation of J/tp photoproduction at HERA has been carried out using values of the 
nonperturbative colour-octet terms determined from a fit @] to the CDF data . These 
calculations predict a cross section for HERA rising with z, which is not seen in the data. 
This is illustrated in Fig.^ where the dashed line shows a sum of the colour- singlet and 
colour-octet contributions both calculated at leading order. 

■^The NLO and z distributions for the ip' have large theoretical uncertainties and cannot be ac- 
counted for simply by a scale factor. 
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9 Conclusion 



We have measured inelastic J/t/) photoproduction in the range 50 < W < 180 GeV and 
0.4 < z < 0.9. The cross section rises with W. In this z interval the photon-gluon fusion 
process is expected to dominate. A NLO calculation for the photon-gluon fusion process 
agrees with the data both for the integrated cross section and the differential distributions 
over z and p'^ in the kinematic range z < 0.8 and p"^ > 1 GeV^, using gluon distribution 
parametrizations compatible with those determined from the F2 measurements performed 
at HERA. The predictions of a specific leading order colour-octet model, as formulated 
to fit the CDF data on J/'^ hadroproduction, are not consistent with the data. 
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W range (GeV) 


50-90 


90-120 


120-150 


150-180 


0.9 < z < 1 


14±4 


29±6 


27±6 


12±5 


0A< z< 0.9 


67±9 


53±8 


35±7 


26±7 


z < 0.4 


19±6 


J/V' — > e+e 


W range (GeV) 


50-90 


90-120 


120-150 


150-180 


0.9 < ^ < 1 




22 ± 6 


0.5 < z < 0.9 




20±6 


33±7 7±3 



Table 1: Number of events in the various z and W ranges for J/ip /i+ii and e+e . 









3/r ^ 


< 0.9 






w 


< > 


A 


o-ep(0.4< z <0.9) 


aep{z <0.9) 




cr^p(2; <0.9) 


(GeV) 


(GeV) 




(nb) 


(nb) 




(nb) 


50-90 


73 


18% 


1.74±0.23l^:i^ 


1.92±0.25;^:^^ 


0.0555 


34.6±4.5;^:'^ 


90-120 


105 


24% 


i.io±o.i7l^:i^ 


1.23±0.19l^:i^ 


0.0232 


53.0±8.2ll^ 


120-150 


134 


24% 


o.75±o.i4j:^:^^ 


o.84±o.i6j:^:^^ 


0.0157 


53.5±io.2j:^:^ 


150-180 


162 


16% 


o.8i±o.2oj:^:^? 


0.90±0.22l^:i^ 


0.0110 


81.8±20.0l^:^ 










-, z < 0.9 






W 


<W> 


A 


(Tep(0.5< z <0.9) 


aepiz <0.9) 




a^p{z <0.9) 


(GeV) 


(GeV) 




(nb) 


(nb) 




(nb) 


90-120 


107 


16% 


0.63±0.20l^:i^ 


0.79±0.25l^:it 


0.0232 


34.i±io.8j:;^:^ 


120-150 


136 


20% 


0.85±0.18l^:i^ 


1.05±0.22l^:i^ 


0.0157 


66.9±14.0l^J 


150-180 


166 


7% 


0.55±0.25l^:^^ 


o.68±o.3ij:^:^^ 


0.0110 


61.8±28.2l^5^o 






J 


/r it+ij-, z < 0.8, p\. > 1 GeV^ 






W 


<W> 


A 


(7ep(0.4< z <0.8) 


aepiz <0.8) 


$r 


a^p{z <0.8) 


(GeV) 


(GeV) 




(nb) 


(nb) 




(nb) 


50-80 


70 


16% 


0.80±0.18l^:^^ 


0.93±0.21+^:i^ 


0.0452 


20.6±4.6l^:^ 


80-110 


96 


22% 


0.65±0.12l^:^^ 


0.76±0.14l^:^? 


0.0268 


28.4±5.2+^;^ 


110-180 


138 


18% 


0.89±0.19l^:^^ 


1.04±0.22l^:i^ 


0.0334 


31.1±6.6l^:5 


J/^ ^ e+e", ^ < 0.8, > 1 GeV^ 


W 


<W> 


A 


C7e3,(0.5< z <0.8) 


aep{z <0.8) 




ajp{z <0.8) 


(GeV) 


(GeV) 




(nb) 


(nb) 




(nb) 


110-180 


142 


15% 


0.70±0.18in9 


0.96±0.25l^i^ 


0.0334 


28.7±7.5l^:^ 



Table 2: Inelastic J/tp photoproduction cross sections for the muon and electron decay 
modes in the two regions z < 0.9 and z < 0.8, > 1 GeV^. From left to right we give 

the W range, the W mean value {< W >), the acceptance (A), the measured ep cross 
sections, the ep cross sections extrapolated to ;z = 0, the flux factor ($t) and finally the 
7P cross sections. 
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w 

(GeV) 


<W > 
(GeV) 


a^p (0.4 < z< 0.9) 

(nb) 


a^p {z < 0.9) 
(nb) 


50-90 


73 


31.4±4.ll^:t 


34.6±4.5lt:i 


90-120 


106 


41.4±6.0l^:^ 


46.1±7.0l^;^ 


120-150 


135 


5i.o±7.6j:t:9 


57.8±8.7l^:^ 


150-180 


162 


68.2±14.5l^:| 


75.3±16.6j:^:^ 


W 
(GeV) 


<W > 
(GeV) 


a^p (0.4 <z< 0.8, p| > 1 GeV^) 
(nb) 


cr^p {z < 0.8, p| > 1 GeV^) 
(nb) 


50-80 


70 


17.7±4.0+^:^ 


20.6±4.6+^i 


80-110 


96 


24.3±4.5l^:^ 


28.4±5.2l^;^ 


110-180 


140 


25.7±4.2+^:t 


30.1±5.2t^:^ 



Table 3: Cross sections o^^^jj^x for the phase space regions: 0.4 < z < 0.9 (top-left), 
z < 0.9 (top-right), 0.4 < z < 0.8 and > I GeV^ (bottom-left), z < 0.8 and Px > ^ 
GeV^ (bottom-right). The cross sections for 0.4 < z < 0.9 and for z < 0.9 in the W range 
from 50 to 90 GeV come from the muon channel only. The first error is statistical, the 
second one comes from all the systematic errors added in quadrature. The other three 
measurements come from the combination of the electron and muon results as described 
in the text. The first error contains the contribution from statistical and decay channel 
specific errors while the second contains all sources of common systematic errors. In the 
regions 0.4 < 2; < 0.8, > 1 GeV^ and z < 0.8, > 1 GeV^ for W in the range from 50 
to 110 GeV only the muon channel data are used while for the highest bin electron and 
muon results were combined as explained in the text. 
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Figure 1: Dominant inelastic J/^ production mechanisms at HERA. Photon-gluon fu- 
sion is described by diagram (a). Diagrams (b) and (c) correspond to diffractive proton 
dissociation and resolved photon J/'^ production, respectively. 
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Figure 2: Invariant mass spectrum for the muon pair sample (a) for 0.4 < z < 0.9 and (b) 
for z < 0.4, in the W range 50 to 180 GeV. The invariant mass spectrum for the electron 
pair sample (0.5 < z < 0.9 and 90 < 1^ < 180 GeV) is shown in (c). The muon mass 
spectrum (a) was fitted to the sum of a Gaussian and a flat background; the spectrum (b) 
was fitted to the sum of a Gaussian and a linear background. The electron mass spectrum 
(c) was fitted to the sum of the convolution of a Gaussian and a bremsstrahlung function 
plus a linear background. 
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Figure 3: In (a) the uncorrected W distribution of the /i+yU" data (full dots) with 0.4 < 
2; < 1 is compared to the mixture of HERWIG and EPSOFT (continuous histogram) 
described in section ^ In (b), (c) and (d) similar comparisons between data and the 
Monte Carlo mixture are shown for the distributions of px and polar angle 6 of the J/ip 
and for the energy in the forward calorimeter EpcAL, respectively. The Monte Carlo 
mixture is normalized to the number of measured events. 
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• ZEUS 
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O FTPS 
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MRSG (m,= 1.4 GeV, AocD=300 MeV) 

CTEQ3M (m,= 1.4 GeV, Aocd=300 MeV) 

GRV (m,= 1.55 GeV, Aoco=215 MeV) 

z < 0.9 



50 



100 



150 



200 

W (GeV) 



Figure 4: The direct inelastic J ftp photoproduction cross section as a function of W for 
Data from ZEUS, HI 



z < 0.9. Data from ZEUS, HI [y], FTPS |[, NA14 |Tg and EMC are shown. 
The ZEUS result at the lowest W value is obtained with the muon channel only. The 
inner error bar indicates the statistical uncertainty, the outer error bar the quadratic 
sum of the statistical and systematic uncertainties. This is also true for the results from 
HI, FTPS and EMC Collaborations. The other three ZEUS measurements come from 
the combination of the electron and muon results as described in the text. The inner 
error bars represent the statistical and decay channel specific errors added in quadrature, 
the outer ones the statistical, decay specific and common systematic errors added in 
quadrature. The lines correspond to the NLO prediction from [Q] assuming the GRV p7 
(continuous), MRSG (dashed) and CTEQ3M Q (dotted-dashed) gluon distributions 
with nic = 1.4 GeV and Aqcd = 300 MeV, the dotted curve was obtained with GRV, 
rric = 1.55 GeV and Aqcd = 215 MeV . The curves are scaled up by a factor of 1.15 to 
take into account the contribution from tp' — > J/tpX. 
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Figure 5: The direct inelastic ^/ip photoproduction cross section as a function of W 
for z < 0.8 and > 1 GeV^. Data from ZEUS and HI ^ are shown. The ZEUS 
results in the lowest two W bins are obtained with the muon channel only. The inner 
error bars indicate the statistical uncertainties, the outer error bars the quadratic sum of 
the statistical and systematic uncertainties. This is also true for the results from the HI 
Collaboration. The ZEUS measurement in the highest W bin comes from the combination 
of the electron and muon results as described in the text. The inner error bar represents 
the statistical and decay channel specific errors added in quadrature, the outer one the 
statistical, decay specific and common systematic errors added in quadrature. The lines 
correspond to the NLO prediction from assuming the GRV (continuous), MRSG 
Q (dashed) and CTEQ3M H (dotted-dashed) gluon distributions with = 1.4 GeV 
and Aqcd = 300 MeV. The curves are scaled up by a factor of 1.15 to take into account 
the contribution from ip' —>■ J/ipX. 
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Figure 6: Differential cross section da/dp\ for tfie inelastic ^/ip — > /i^yU^ sample with 
50 < < 180 GeV and z < 0.9. Data from ZEUS and HI ||15| are shown. The error 
bars indicate the quadratic sum of the statistical and systematic uncertainties. The NLO 
computation with the GRV structure function, rric = 1.4 GeV and Aqcd = 300 
MeV is shown as the solid line. The theoretical curve is drawn only for > 1 GeV^ 
because in the low px region the calculation is not reliable. In the theoretical curve the 
15% contribution of the ip' has not been included. 
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Figure 7: Differential cross section da/dz for the inelastic ^/ip — 
50 < < 180 GeV and > 1 GeV^. Data from ZEUS and HI 
inner error bars indicate the statistical uncertainties, the outer error bars the quadratic 
sum of the statistical and systematic uncertainties. The NLO computation |T| with the 
GRV structure function, rric = 1.4 GeV and Aqcd = 300 MeV is shown as a solid 
line. The dashed line is given by the sum of the colour- singlet and the colour-octet leading 
order calculations In the theoretical curves the 15% contribution of the ip' has not 
been included. 
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